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Abstract 



The self-consistent statistical approach to the problem of planetary and 
stellar magnetism is suggested. The mechanism of magnetic fields generation 
in the astronomical objects, where the existence of fields is associated with 
the axial rotation of objects, is discussed. In the general case the light 
pressure, the centrifugal, gravitational and other forces produce partial p- 
separation of the charges. As a result of the system rotation, the magnetic 
00 ■ fields of the currents of these charges are not compensated. The influence 

^> ' of various factors on the magnetic field of some object is analysed. 

o ' 
o 



1 Introduction 

The origin of magnetism in various astronomical objects is of great theoretical 
and practical interest. The existence of magnetic fields in these cases indicates 
that currents are circulating inside such objects. What is the nature of these 
currents? 

At the present time there is a great deal of information on this subject. In 
^ ■ the literature the existence of magnetic fields is associated with: 

1) various dynamo processes: a) the classical kinematic dynamo-models 
(Busse, 1979; Parker, 1979; Moffat, 1978); b) the kink instability mechanism 
(Alfven et al., 1974; Alfven, 1981); c) the action of gravitational forces (Hide, 
1956); d) the specific models (heat convection, precession) (Urey, 1952; Busse, 
1976; Bullard et al., 1971; Runcorn, 1975); 

2) or with the remanent magnetization (Sharpe et al., 1976). 
These theories encounter some problems: (Runcorn, 1975; Malkus, 1963; 

Stevenson, 1974); the Earth's core paradox (Higgins et al., 1971; Kennedy et al., 
1973); the energetic contradictions (Jacobs et al., 1972); the contradictions of 
precession mechanism (Rochester et al., 1975; Busse, 1971) etc. 

A variety of approaches to the problem of magnetism of astronomical objects 
may exist. One of them consists in describing the process of magnetic field birth 
and generation up to observable values. This approach is undoubtedly important 
and interesting, from principal positions, to elucidate the mechanism of such a 
generation. As an example we refer to the "battery theory" for the generation 
of cosmic magnetic fields (Biermann 1950; Mestel et al., 1962; Dolginov, 1988). 
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Note, that the paper is concerned with the above mentioned theory (in initial 
cause). However, some differences exist. The basis for the "battery theory" 
is the "Biermann effect", where the following term provides the magnetic field 
generation: 

c VP e 1 
cVr/ x VT = --V x -^4 = -M,V x g? fJ , 

where r\ is the thermopower constant, P e is the electron pressure, g^ff = 
g + ilxilxris the sum of gravitational and centrifugal accelerations (the 
causes of generation: nonuniform rotation, the inhomogeneous distribution of 
chemical elements, temperature inhomogeneities) . At the same time, however, 
it has rather essential disadvantages. First, one must proceed from some un- 
known hypothetic state. Second, for the time during which the magnetic field 
is increasing, so many various random events take place, that they can not be 
taken into consideration even in principle. For example, random flows of charged 
particles (cosmic rays) can violate the 2™ magnification (as a consequence of 
long-range Coulomb interactions). Third, even the main factors can be taken 
into account within the framework of approximate models only (the continuum 
medium model, for instance). Besides, there exist problems for calculation. At 
present it is impossible to solve the full system of equations (Dolginov, 1988). 
The magnitude of the stationary field cannot be determined from the linearized 
equations (note that the nonlinear terms are the same order of magnitude as 
the linear ones). But investigation of nonlinear effects are still in their infancy. 
All these difficulties do not allow quantitative characteristics of the field to be 
determined. There exists, however, some another approach to this problem. It 
is as follows. At present, a set of observable characteristics exists for studying 
astronomical objects. These characteristics (seem to be stable) include such ones, 
as the frequency of rotation of a system, its composition, electric and magnetic 
fields, etc. Here, the question arises: is there any connection between these ob- 
servable quantities, so that, proceeding from some particular observable system's 
parameters, one can calculate other parameters. 

Many astronomical objects (nebulae, planets, the Sun and other stars) are 
known to be rotating bodies. The attempt to find the relation between these 
phenomena is made in this article. Since some astronomical magnetic fields are 
rather stable, it can be supposed stable or quasistable state exists in the system 
and the microscopic statistical approach to this problem can be proposed. Devel- 
oping this approach, an attempt is made in this paper to describe theoretically 
stable states of the system and the dynamics of fields. In principle, the statis- 
tical approach allows the exact solution for the stationary field to be obtained 
in integral form. Some distinct value of the field (to sufficient accuracy) can 
be calculated from this solution for all specific values of parameters ( nonlinear 
effects are implicitly taken into consideration). 

The peculiarities of this work are as follows: 

- the mechanism dispenses with the need for initial magnetic fields; 



2 



- the mechanism proposed can explain the existence of magnetic fields in a 
wide majority of astronomical objects, i.e. it is rather universal; 

- the mechanism has self-consistent characteristics and leans upon the micro- 
scopic properties of matter, rather than upon the continuous medium models; 

- the mechanism is simple enough, that is, there is no need to make any 
additional assumption on properties of the inner region of objects (good conduc- 
tion, for example) or on processes occurring inside the objects (convection, for 
example) . 

Note, that the suggestion of the statistical approach to the problem is the most 
important part of the article, since there exist a possibility to study magnetic 
fields without knowledge about the nature of electric double layers: the statistical 
method can include different mechanisms which can lead to separation of charges 
(but not only suggested). For example, the reverse problem can be formulated: 
knowing fields, to find the possible position of electric double layers and theirs 
magnitudes. After that the discussion of magnetism (concrete mechanism) is 
more objective. However, it is not our intention to follow this way in the article. 

In section 2.1 the separation of charges and the magnetic field appearance are 
investigated for the rotating plasma model. This model can be applied to some 
nebulae, for example. In section 2.2 the action of a gravitational force in the 
compact rotating system is additionally considered. This model can explain the 
existence of stellar magnetism. In section 2.3 the magnetism of planets (some 
"cooled" objects) is discussed. In this case the "cooled" (magnetic) materials 
in the surface layer of a planet (the Earth crust, for example) at temperature 
less than the Curie temperature can influence the magnetic field. In section 2.4 
the general remarks on the Kauling theorem and different mechanisms, and some 
evaluations are presented. Section 3.1 outlines self-sustained rotation of a system 
and the equilibrium configuration of electric and magnetic fields. Section 3.2 
analyses the influence of various internal and external factors and the dynamics 
of fields. 



2.1 The magnetic field in rotating plasma 

The mechanism of magnetic field generation in the astronomical objects can be 
qualitatively explained in terms of this model. The rotating plasma model can be 
applied for explanation of magnetic fields existence in rotating plasmoids, such 
as ionized cloud, hot nebulae, etc. 

As a consequence of long-range interaction between the charged particles, the 
plasma essentially differs from gas in some respects. The rotating plasma has 
some peculiarities as compared to rotating objects involving neutral particles. 
The existence of the centrifugal force 
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pe P , 
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(where fi is the angular frequency of the system, m is the particle mass, p is 
its distance from the axis of rotation) causes different effects on particles of 
different masses. As a result of these different effects, the p-dependences of 
particle concentrations are bound to be different for particles of different masses. 
According to the Boltzmann distribution, 

n a = n 0a exp(-^), (2) 

where T is the system temperature, k is the Boltzmann constant, no a is the 
particles concentration on the axis of rotation, U a is the potential energy of 
particles of a-sort; for neutral particles 

U. = -=S^. (3) 

Since the particles of different masses in plasma have different charges, this par- 
tial separation of particles produces a partial separation of charges, that is, the 
negatively charged region must lie near the axis of rotation, whereas the posi- 
tiveLy charged region must lie near the system boundary (the distance Ro from 
the axis of rotation). The electric field Eo(p, z) (in the polar coordinate system) 
exists as a result of this separation of charges. 

This field opposes the considerable separation of charges. Therefore, the 
effect of partial separation of charges is bound to be often negligible for small 
systems, whereas the separation of charges can exert considerable influence on 
some physical characteristics of astronomical objects. 

There is some distance from the axis of rotation Ri(z), where the local charge 
equals zero (some transverse section to the axis). For < p < R\{z) the plasma 
is negatively charged on the average; for R\(z) < p < Ro the plasma is positively 
charged. Therefore, the field Eo(p, z) at i?i(0) is directed towards the axis of 
rotation. 

Of mine interest here is the fact, that the partial separation of charges gives 
rise to the p-dependence of charge density, that is, the plasma possessing a given 
charge density moves round a circle of distinct radius. Therefore, the circulating 
currents are inside the rotating plasma despite the fact, that all particles revolve 
with fi on the average (fi has a distinct direction). For R\ < p < Rq the 
current flows in the direction of rotation (since the charge is positive), whereas 
for < p < Ri the current is opposite to the direction of plasma rotation (since 
the charge is negative) . In the general case the magnetic actions of these currents 
are not compensated and the magnetic field exists. This is the magnetic field of 
a solenoid whose axis coincides with the axis of rotation of the system. 

It follows from the Biot-Savart law, that the magnetic field is 

H i U' xR '' 



c J(y ) R 6 
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where V is the volume of the system, c is the speed of light, dv' is the region 
with currents, R' is the distance of this region to the point of observation, j' is 
the current density. We have 

R' = (R p ; Rp, R z ) = (p - pi cos((^i - ip); pi sin(^ -<pi);z- zi); 

dv' = p 1 dp 1 dtp 1 dz 1 , j' = (0; , j<p(pi, <pi, z±) = etlpiq(pi, <pi, Zi). 

The electric charge density q is 

q(Pi,zi) = ^exp (^--j^f u ij ~ «o e exp ^~^ U ej , (5) 

For rotating, fully ionized plasma, which consists of the elements with atomic 
number TVj, it follows in polar coordinates p,(p,z (axis z is in the fi direction), 
that 



jj = _mn 2 p( 



2 



■>2„2 

Ui = - 



+ e f Pl E 0p (p', Zl )dp' + e / 1 E 0z (pi,z')dz'+ 
Jo Jo 

rpi P ff rzi p ff 

/ e Jl±VLp'dp' - / ^n Pl dz' , (6) 
Jo c Jo c 

M£i z ni I'f 1 f zi 

— — ^ - eiV, / E 0p (p', Zl )dp' - eNi / E 0z ( Pl , z')dz'- 
£ Jo Jo 



/ ^sip'dp' + / ^n Pl dz' , (7) 

Jo c Jo c 

where n e (p,z) is the concentration of electrons, rii(p,z) is the concentration of 
ions, — e is the electron charge (e > 0), m is the electron mass, M is the ion 
mass, the electron and ion temperatures T e = T% = T; the quantities riQ a can be 
found from conditions: f,y\ n a dv' = N a K a , where K a is the total number of 
particles of a-sort, N e = 1 and K e = NiKi + A, A is the electron surplus. 
The electric field can be found from the integral equation 



Jv 



eR'(p',z')q(p',z')dv' 



(v) R' 3 

Heat movements counteract the compensation of the charge density: the 
statistical (dynamical) equilibrium is established. In the general case for ro- 
tating plasma of complex composition, — > Y^i3% j n 0a — * n^ a . In prin- 
ciple, we can introduce the "effective" ion mass M^^ (in the general case 
M- ff > Mf v = (J2MiKi)/(J2Ki)) and N- ff . Note, that a small quantity 
of heavy ions plays an important role. As follows from the derived expressions, 
the effect is proportional to the particle's concentration (or, the electron sur- 
plus A). Since for Q = and for Q — > oo the effect vanishes, there exists such 
a value Cl m that | H z ($l m ) \= max. Analogously, there exists T m for which 
I EL z (T m ^ |= Txidx. The dependence H z (yi,T^j can be represented as a "pit". This 
is only the discussion of principles of approaches to the magnetic field generation 
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problem. Strictly speaking, the model uses the externally given fixed "wall". To 
use the described model for real objects, the following factors need to be taken 
into consideration: the gravitational force, if it is important (maintains the sys- 
tem as a whole, for example), or, the expansion of system boundary (the change 
of integration limits, if the equilibrium distribution has a chance to be set) and 
the dependence tt(p,<p,z) (nonrigid rotation). Some practical remarks will be 
given in Section 2.4. 



2.2 The magnetic fields of rotating stars 

The consideration of the starry magnetism differs from the case discussed above 
both quantitatively (ft and no a may be considerably larger than corresponding 
quantities for the ionized cloud) and qualitatively. One of such qualitative differ- 
ences lies in the fact, that the gravitational force must be taken into account (the 
external boundary doesn't be fixed, but can be determined in the self-consistent 
manner) : 

„ M(r)mr , , 

F g = , (8) 

where 7 is the gravitational constant, r is the distance from the center of a star, 
M(r) is the mass of a star part in the < r' < r region of the star (the star is 
spherically symmetrical) . 

The gravitational force is directed toward the star center, that is, this force 
has the /5-projection (toward the axis of rotation). Therefore, this force competes 
with the centrifugal force (the actions of these forces are opposite) in the processes 
of separation of charges and generation of electric and magnetic fields. It is 
clear from the general reasoning, that at the given angular frequency and starry 
composition the magnetic field reverses its direction at some "critical" mass (or at 
"critical" Rq c ); or the magnetic field of a certain star reverses its direction at some 
certain critical angular frequency £l c . In the case of gravitational confinement, 
the (more light) electron cloud is distributed near the boundary (and just near 
the axis) , but the (more heavy) ion layer is distributed near the center (the value 
kT is the same for all particles and « U g ). 

The potential energy of some particle with mass m is: 

U 9 = , (9) 



where r = \J p 2 + z 2 . 

In the general case the distribution of charges includes three charged regions. 
The negatively charged region is placed near the axis of rotation and the system 
boundary; the positively charged region (torus) is placed near (and far) the radius 
which can be determined from the condition F g = F c . The bounds between the 
charged regions can be found from the equation n e (p,z) = J2i n i{p, z )- F° r 
example, it follows that the equation for plasma consisting of particles with one 
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number N is: 

no e exp 



kT 



f z eH [P eH f 

/ — f -n Pl d Zl - / — -fWpi - e / E' dr'+ 
Jo c Jo c J 



7" 



n 0i exp [ 



kT 



m e Q p M(r)m e 
2~ 

MiQ 2 p 2 



+ Ne 



eH r 



Qpidzi — 7 



/ E '°*' + i 

M(r)M, 



{Ipidpi — 



(10) 



In the general case the components Eo(p, z) can be found from the system of 
integral equations: 



Eo P {p,z) = e 



(V) 



pi[p - picos(ipi - ip)]dp 1 dz 1 dip 1 



Eoz(p,z) = e 



q(pi,zi) 



[z — z\)p\dp\dz\dtp\ 



i(v) R' 3 

R' = \J(z- Z]) 2 + p 2 + p\- 2pp 1 cos(ip 1 - ip) 

with additional term jM(r)Mi/r in Ui and ^M(r)m e /r in U e . 

The magnetic field can be obtained from the following expressions: 



(11) 
(12) 



1 



3<p{pi,zi) 



H z = - 

cJ(v) 



c J (V ) R' 3 
3tp(pi,zi 



(z — z\)p\dp\dz\dtp\ , 



R>3 



pi[picos((pi - ip) - p]dpidzid(pi 



(13) 
(14) 



In the case of unsteady rotation (Q 2 p 2 /2) — > f£ft(p', z)p'dp' . 

As a first approximation, the effect is proportional to the particle's con- 
centration (or, the ion surplus —A). There exists such a value T m , that 
I H z (T m ) |= max. At first, the value H z increases (to a maximum iT^f^J > 0) 
with increasing the value Q to furthermore, the value H z decreases (to a 
minimum H z (Q 2 n ) < 0) with increasing the value 0, and, finally, the value H z 
tends to zero with further increasing Q. We can write approximately 



dH z H z 



n dH z 



dT 2T 2T dn ' 

The dependence H Z (Q,T) can be represented as "a hill passing to a pit". Anal- 
ogously, the dependence H Z (M,T) can be described as "a pit passing to a hill". 
However, in the real case the value T can increase with increasing the body mass 
M, and the value H z can be finite (^ 0) with increasing the value M. 
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There is another qualitative difference between the cases of star and rotating 
plasma. The star is a nonequilibrium system. At some point of the star the 
radiatmon from and to the center of the star are not compensated, thus causing 
the light pressure other than zero to exist. The action of this light pressure 
for particles depends on the effective section of interaction of these particles. 
Therefore, tie light pressure force for ions is larger, than that for electrons, and 
this influence on the charge separation process is opposite to the gravitational 
force action. 

Given the effective sections of interactions a e and oi and the light pressure 
function P(p, z), the field H can be obtained in terms of this effect. To do this, 
the following substitutions need to be done: 

-m e n 2 p 2 -m e Q 2 p 2 f(fi' z ) 



2 2 7(o,o) 

-M$l 2 p 2 -M,iVi 2 p 2 fM 



[ {p ' z> i i i 
- o e / P(p ,z)dr; 

J (0,0) 



2 2 7(o.o) 

Some important remarks will be made in section 2.4 



a t / P(p',z')dr'. 
J(o.o) 



2.3 The planetary magnetism 

In this case the radiation can be ignored. In some instances the gravitational 
force can be ignored only when the of charge separation region is small, or, 
more precisely, its action is reduced to the action of pressure on the boundary 
of this region, ft and iVj quantities for the planets may be larger, than those 
for the other astronomical objects. The distribution of electron (or ion surplus) 
is determined by characteristics of planetary system (Boltzmann factor depends 
on the solar mass, chemical composition etc). For some planets with strong 
magnetic field the radiation belts (the Earth radiation belt, for example, which 
is the rotating plasma substantially) have influence on the magnetic field. Note, 
that radiation belts can be considered as a separate system (with their own 
boundary conditions), for which all derived formulae are applicable. 

The case of planets has some qualitative difference from the other cases. The 
planets are relatively cold objects. Therefore, two remarks are in order. First, 
not all particles are involved in the processes of charge separation and of the 
magnetic field appearance. The neutral particles influence on the electric field 
Eo(p, z). To take into account this influence, the following substitution needs to 
be done in the denominators of equations (11), (12), (13), (14): R' 3 — > R' 3 e(p,z), 
where e(p, z) is the permittivity. The polarization of nonmagnetic materials can 
influence on the magnetic field. 

Second, the influence of magnetic materials in the surface layer of a planet 
on the magnetic field at temperature lower than the Curie temperature must be 
taken into account: 

B = pH. + 4vrM , 
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where H is the magnetic field from (13) or (14), p, is the permeability, Mo is the 
remanent magnetization. 

The latter factor can become the only one which maintains the existence of 
the magnetic field as the axial rotation of a planet is decelerated. The qualitative 
picture, described above in previous subsection, is valid for the charge distribu- 
tion and the dependence of magnetic field on different factors. Note, that (an 
agreement with this paper) there exists the electron surplus near the Earth's sur- 
face (see The Pysical Values, 1991). The most significant remarks with practical 
viewpoint is made in the next subsection. 

2.4 Some remarks and estimations 

First of all we note that the Kauling theorem does not valid for models described. 
The electric field possesses the E$ p and Eq z components only (Eq v = 0). However, 
the nature of currents does not determined by the magnetic force ^[v x B]. 
The existence of currents is the consequence of the system rotation: portions 
of charged particles (which form charged regions) are involved in the rotation. 
In spite of particles movement (rotation) with the same ft (without friction), 
currents can have the mutually reverse directions. There is a surface with j = 0; 
the field B can be equal zero on some surface; rot B can also be equal zero. 
Thus, the principal conclusion consists in the fact, that the magnetic field can 
origin in a system with cylindrical symmetry. 

Second remark needs to be made in connection with conditions of application 
of models. It is easy to make an estimate of the field produced by the effect 
described for T = 0. The magnetic field (in Gaussian units) is H ~ RJ/c, here R 
is the system radius, J ~ pv. If the gravitational force is most important, then 
eE ~ Mj</, p ~ Mj</ / (eR) here g' ~ jpbR, where j is the gravitational constant, 
pi, is the body density. Substituting typical numbers (Mj ~ 10 _23 gram, v/c ~ 
lCT 4 ,e ~ 10~ 10 esu, p\ } ~ 10 _3 gram/cm 3 ), we have H ~ 10~ 4 Gauss for R ~ 
10 9 a.u. However, the model is statistical (principally): statistical (dynamical) 
equilibrium holds by heat movements, and the originated electric force cannot 
compensate the charge density (for example, the noncompensated force F g exists, 
but the earth atmosphere doesn't "fall" on the Earth). The estimations needs to 
be corrected for stellar and planetary magnetism. To do this, two ways exist. 1) 
The above estimations suppose that the magnetic field is produced by a charge 
surplus. Note, that if for an earthlike body H ~ 1 — lCT 4 Gauss, then the surplus 
charge density p needs to be equal p ~ 10~ 4 — 10 _8 esu/cm 3 , that is the number 
of charged particles is 10 2 — 10 6 per cm 3 . It is possible for the earth phenomena 
(see The Physical Values, 1991). 2) The "elucidation" of the gravitational force 
for system neutral as a whole does not meant that the body will scatter under 
influence of electric force. These forces do not compensate. Besides, any liquid 
or rigid body can possess rather great portion of charges without interruption. 
The "capacitor" (described above the double layer) discharges. To account the 
latter effect, one can use some potential (barrier maintaining a "contact potential 
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difference") which counteracts the discharge process. Note as a remark, that 
Earth structure (in chemical and material composition) can be determined by 
self-consistent manner (using the chemical potential fi). The presence of charged 
particles can be described as a "dilute solution", and the equilibrium condition 
for different phases is: p, + U = const. 

The following remark concerns the role of electromagnetic forces for Earth 
phenomena. As a consequence of long-range (rather strong) electric interaction, 
some small charge surplus is suffice for comparability the electric energy with the 
motion energy in atmospheric (different winds) and hydrodynamical (sea cur- 
rents) phenomena. To account this, one must additionally write: 1) the motion 
equation for charged particles, and, 2) the equation for the interaction of charged 
particles and neutral gas (or liquid) masses. Cosmic rays influence on the earth 
activity (earthquake, vulcan activity), since the presence of magnetic materials 
(hysteresis) leads to heat production in inner earth regions. 

3 Magnetic fields and some factors 
3.1 Self-sustained equilibrium rotation 

It is shown above, that in a system possessing some fraction of free charges a par- 
tial separation of charges occurs under an effect of centrifugal and gravitational 
forces (n e (p, z) does not coincide with rii(p, z)) and, as a result, the electric field 
E(p, z) arises. Due to system's rotation the magnetic effect of currents of these 
charged regions is not mutually compensated, and the magnetic field H arises. 
In this case the rotation frequency Q(p, z) was assumed to be specified. 

The statement of the problem may be changed. For this purpose we notice 
that the field H is directed along the axis of rotation (axis z), and the field E has 
two components E p and E z only, i.e. it has a component directed towards the 
system's axis of rotation. As a result, the drift of particles will take place in cross 
fields Eo p e p and He z . This drift motion will occur over circles around the axis 
of rotation. Hence, it is worthwhile to find such a rotation rate, that resulting 
electric and magnetic fields lead to a drift rotation with the same rate, i.e. one 
should find stable rotational states of a charged particles in a self-consistent 
manner. Note, that the (averaged) movement of charged particles must be close 
to the rotation of the system as a whole (a consequence of the system stability). 

From condition = Qp and the expression for drift velocity 



e<pcE 0p (p,z) 
H(p,z) 



the unknown rate is determined as 



n = 



cE 0p 



(15) 



p[H (p, z) + H i y 
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where p is the distance to the axis, c is the light velocity, H\ is an outer or 
intrinsic magnetic field (Hi > Hq), and quantities Eq p and Hq can be taken from 
previous section. 



3.2 Dynamics of magnetic fields 

Let us imagine that the external magnetic field Hi suddenly increases: Hi — > 
Hi + H'i. As a result, as follows from (15), the equilibrium frequency of rotation 
of the object (or a fraction of charged particles) decreases in this case. This 
decrease of £l(p, z) results in decreasing the field Ho — > Ho — H and, hence, the 
resulting field slightly lowers: H = Ho + Hi + H' x — H . This is a consequence 
of the general Le Shatelier principle. Similarly, the weakening of the external 
magnetic field Hi leads to increasing the drift rotation rate, which, in its turn, 
causes the increase of the field Ho, i.e. the processes weakening the external 
action occur again (the system tends to conserve its equilibrium state). 
Note that at equilibrium either $7 = const, or Q = £l(p, z): 

Q 2 n 2 CP 

^ - / n(p',z) P 'dp', (16) 

l Jo 

for charged particles and the source of energy and the mechanism of compensation 
of heat losses (for friction) must exist. If the state of a system differs from 
equilibrium one by any reason (for example, the external forces are acting, or 
some processes occur inside the object, which lead to redistribution of charges), 
then, in principle, one can find the forces, tending to bring the system into 
equilibrium, and describe the dynamics (of fields, of rotation rate). 

It follows from (15), that as the rotation frequency of an astronomical object 
changes adiabatically, the electric and magnetic fields also change: 

o _ c \ E op _ E 0p^_ ] 
p H H 2 J ' 

As the first approximation, using 

E = | 9E °p^ | dE Qpf, 9E °r> - 2TdE °? 



on OA &T ' on n or 



we have 




H = H 1^-7; T-2T- . (17) 



The value ^jf- must be evaluated from the self-consistent system of equations 
(from Sec. 2). Since H = Hq + Hi, the hysteresis influences on the field dynamics. 
The change of the magnetic field polarity may occur (without changing the direc- 
tion of astronomical object's rotation), when the direction of the electric field Eq p 
inside the system changes. The change of polarity for a system with magnetic 
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materials (there exists some threshold due to hysteresis) requires long-lasting 
influence (cosmic rays can change the total charge of the system, for example). 

In principle, one can obtain the solution to the problem of fields dynamics, in 
the case of "sufficiently slow" variation of various external and internal factors. 
It is apparent that the dynamics of T, p, (V) and Q determines the dynamics of 
electric and magnetic fields by formulae (15) - (17). Therefore, to finally deter- 
mine the fields dynamics and the rotation rate, one should also add the following 
equations: 

1. the moment of external forces acting on a system determines the additional 
change of rotation frequency £1 (if the object does not possess central symmetry); 
this factor may cause the magnetic field precession; 

2. the law of variation of object's form (V) (this factor is insignificant for stars 
and planets); 

3. the law of variation of density p(R, t); 

4. the law determining the temperature variation T(R, t); 

5. the law of (external) change of the total charge (this factor is the most im- 
portant one). 

Now it is clear that the characteristic times of variation of all these quantities 
must be higher than characteristic times of diffusion for a fraction of charged 
particles participating in the fields generation. In the opposite case one must take 
into account diffusion terms in equations for p, T, and (V) (i.e. the equation 
of diffusion should be written separately for a fraction of charged particles). In 
this case the hysteresis phenomena are characteristic of the fields variation. 

How the object driven from the state of equilibrium will behave? From gen- 
eral considerations, the transition into a new equilibrium state (or into in old 
one, if the factor acts for a limited time) may occur either asymptotically, or in 
an oscillational manner (with polarity change, in particular), depending on the 
fact, whether the characteristic time of system oscillation (^ _1 ) is higher (or, 
respectively, considerably lower) than that of system's energy dissipation. 

In the case, if the rotation frequency inside an astronomical object is lower 
than the equilibrium frequency possible for the given system, then the new equi- 
librium state will correspond to the absence of rotation. Because of the de- 
pendence Q(p, z) / const (caused by the fields) the friction will arise between 
different sections of the object and, as a result of energy dissipation, the rotation 
of an astronomical object will be slowed down. This is just an internal mecha- 
nism of slowing down the rotation of astronomical objects (and warming-up of 
inner regions). 

4 Conclusions 

The general statistical approach to the magnetism problem is introduced. It 
can help to obtain the value of the magnetic field to sufficient accuracy (from 
the exact integral solution) . The mechanism of the magnetic field generation in 
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astronomical objects (nebulae, stars, planets) can be explained in the following 
way: in the general case the centrifugal force, the gravitational force and the 
light pressure act on dissimilar particles in a different manner, causing different p- 
distributions of these particles and giving rise to the partial separation of charges; 
the electric field can be found from formulae (11), (12) (see the substitutions in 
the text); the rotation of these objects produces (from these charges) circulating 
currents, the magnetic action of which is not compensated, and the magnetic 
field can be found from formulae (13), (14) (see the substitutions in the text). 
Generally speaking, the existence of the magnetic field dispenses with the need 
for conductive regions; the availability of some portion of free charges would be 
sufficient for this mechanism. 

The equilibrium frequency of rotation of an astronomical object (see (15)) 
may be found in a self-consistent manner. Knowing the factors, which influ- 
ence the density, temperature, rotation of an object and external fields, one can 
determine in this case the dynamics of system's magnetic field. 

The further applications of formulae obtained might be associated with nu- 
merical methods or with some simplifications which take into account some spe- 
cific data. 
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